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The role of charge order in the phase diagram of high temperature cuprate superconductors has been recently 
emphasized by the experimental discovery of an incipient bi-directional charge density wave (CDW) phase in 
the underdoped regime. In a subset of the experiments, the CDW has been found to be accompanied by a d-wave 
intra-unit-cell form factor, indicating modulation of charge density on the oxygen orbitals sandwiched between 
neighboring Cu atoms on the CuO planes (the so-called bond-density wave (BDW) phase). Here we take a 
mean field Qi = (27r/3, 0) and Q 2 = (0, 27r/3) bi-directional BDW phase with a d-wave form factor, which 
closely resembles the experimentally observed charge ordered states in underdoped cuprates, and calculate 
the Fermi surface topology and the resulting quasiparticle Nernst coefficient as a function of temperature and 
doping. We establish that, in the appropriate doping ranges where the low temperature phase (in the absence 
of superconductivity) is a BDW, the Fermi surface consists of electron and hole pockets, resulting in a low 
temperature negative Nernst coefficient as observed in experiments. 


I. INTRODUCTION 

The origin and character of the enigmatic pseudogap phase 
in the underdoped regime of high teniperature cuprate super¬ 
conductors remains an open problenPEI while the insulat¬ 
ing parent compounds of these systems are well understood 
as three-dimensional (3D) antiferromagnetic Mott insulators, 
the normal (non-superconducting) phase above superconduct¬ 
ing transition temperature (Tc) at finite hole doping evinces 
an anisotropic spectral gap (pseudogap) at low energies be¬ 
low a temperature scale T* > Tc and behaves strikingly dif¬ 
ferently from a Fermi liquid. Understanding the pseudogap 
phase from which superconductivity develops at lower tem¬ 
peratures is generally understood to be the key to understand¬ 
ing the d-wave superconducting pairing and the anomalously 
high transition temperature of the superconducting phase of 
the high Tc cuprates. 

Recent theoretical and experimental work has proposed the 
role of various charge, spin, electron nematic, and current or¬ 
dered states competing with superconductivity, and also the 
role of superconducting fluctuations themselves, to explain 
the pseudogap phase above Tc In the class of materi¬ 
als YBa2Cu306+a: (YBCO), exquisite quantum oscillations of 
various electronic properties with the applied magnetic field, 
strong enough to suppress superconducting fluctuations and 
reveal the underlying normal state, have revealed small Fermi 
pockets in excitation spectrum reminiscent of a broken sym¬ 
metry state c omp eting with superconductivity in the under 
doped regime .®ol. The existence of such small Fermi pock¬ 
ets in the underdoped regime combined with a large hole like 
Fermi surface in the overdoped regime, indicates that the nor¬ 
mal state of the cuprates, in the absence of superconductiv¬ 
ity, goes through a Fermi surface reconstruction somewhere 
near optimal doping. Independent evidence of a similar Fermi 
surface reconstruction - from being large and hole-like in the 
overdoped phase to small and electron-like at underdoping - is 
also apparent from the measurements of low temperature Hah 
and Seebeck coefficients which turn from positive at higher 
doping to negative in the underdoped regimJ^^EH, Since the 


signs of the Hah and Seebeck coefficients are determined by 
the sign of the dominant charge carriers, the low temperature 
negative sign of these coefficients in the underdoped regime 
(in the absence of superconductivity) can be explained by the 
existence of electron pockets. Interestingly, the low tempera¬ 
ture Nernst coefficient, which measures the transverse voltage 
induced by a longitudinal thermal gradient in the presence of a 
perpendicular magnetic held, has also been found to be nega¬ 
tive in the under doped regime, while being vanishingly small 
at higher doping. While the sign of the Nernst coefficient, 
unlike that of Hah and Seebeck coefficients, is not directly re¬ 
lated to the sign of the dominant charge carriers, but also de¬ 
pends of the curvature and topology of the Fermi surface, the 
strikingly different behaviors of the low temperature Nernst 
response at low and high hole dopings also point to the exis¬ 
tence of a Fermi surface reconstruction near optimal doping. 
Although various charge, spin, and current ordered states have 
been pro posed to account for the Fermi surface reconstruction 
in ybcoSESI, none had so far been observed in bulk-sensitive 
probes until recently. 

In recent x-ray diffraction experiments two groups have in¬ 
dependently found strong evidence for a short range charge 
density wave phase below the pseudogap temperature scale 
T* for a ra^e of hole doping in the underdoped regime 
of YBCCP^®. In these experiments it is not conclusively 
known if the x-ray diffraction peaks derive from an equal dis¬ 
tribution of domains with uni-directional stripe-like correla¬ 
tions or from correlations with wave vectors (g^i, 0, 0.5) and 
( 0 ,^ 250 . 5 ) (with qi ^ 92 ^ 0.31) co-existing as in a bi¬ 
directional CDW state. However, the lack of anisotropy in 
the scattering signals such as intensities and widths, and also 
salient difference from the stripe like states as observed in the 
LSCO family such as absence of a coincident magnetic or¬ 
der and a strikingly different behavior of the modulation wave 
vector with hole dopin^^, indicate that the short-ranged CDW 
correlations observed in the YBCO family may be different 
from stripes and in fact an incipient CDW order which is bi¬ 
directional. Although the temperature (T) dependence of the 
correlation length above Tc and only short ranged correlations 
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in the Cu 02 planes indicate that the observed charge order is 
only quasi-static, the near divergence of the correlation length 
as T ^ Tc and that the scattering signals significantly in¬ 
crease on application of magnetic field below Tc indicate that 
a true thermodynamic CDW transition at some critical temper¬ 
ature (Tcdw < Tc) may be preempted by the superconduct¬ 
ing transition at Tc. Furthermore, recent inelastic x-ray scat¬ 
tering and nuclear magnetic resonance experiments indicate 
that the sh ort-ra nge charge order observed below T* is in fact 
truly staticP^^, presumably due to pinning by disorder poten¬ 
tial. Evidence for a similar charge density wave transition in 
the underdoped regime has also been found in oth er rep nt 

experiment JSEII. in at least two recent experiment J^lZl the 

charge order has been found to be accompanied by a d-wave 
intra-unit-cell form factor, indicating modulation of charge 
density on the oxygen orbitals sandwiched between neighbor¬ 
ing Cu atoms on the CuO planes (the so-called bond-density 
wave (BDW) state). Taken together, though it is unclear at 
the moment if the bond density wave order observed in the 
cuprates is static and long-ranged or fiuctuating and short- 
ranged below T*, it is clear that its role in the fermiology 
of the cuprates should be significant especially at low tem¬ 
peratures (T < Tc) and in high magnetic fields (sufficient 
to suppress superconductivity) where the bi-directional bond 
density wave is expected to develop long range order resulting 
in Fermi surface pockets in the single particle spectrum. 

An important part of the fermiology of the cuprates is the 
normal state Nernst effect in the pseudogap phase. The Nernst 
response, which measures the transverse voltage induced by a 
longitudinal thermal gradient in the presence of a perpendic¬ 
ular magnetic field, is defined to be positive if dominated by 
vortices in a superconductor. While the quasiparticle Nernst 
signal is typically small for conventional metals due to Sond- 
heimer cancellation, the signal carried by vortices can be large 
and positive in the presence of superconducting fiuctuations, 
as has been found in the cuprates near superconducting Tc 
and above. Suppressing the superconducting fiuctuations by 
a strong magnetic field reveals the normal state Nernst coeffi¬ 
cient (z^/T, with T the temperature) and this has been found 
to drop with decreasing T in the p seudo gap phase, culmi¬ 
nating in a negative z^/T as T ^ (JUMl, xhe low temper¬ 
ature negative Nernst response as T ^ 0 is reminiscent of 
a similar change of sign (with decreasing T) in other trans¬ 
port signatures of the pseudogap phase such as Hall and See- 
beck coefficient The signs of the Hall and Seebeck co¬ 
efficients are determined by the sign of the dominant charge 
carriers and can be explained by the existence of an electron 
pocket centered at ((5/2, (5/2) where the bi-directional BDW 
state is a superposition of CDWs (with d-wave form factors) 
with ordering wave vectors ((5, 0) and (0, Q). This is similar 
to the recently found result of a change of sign (with decreas¬ 
ing temperature) of the Hall and Seebeck coefficients in the 
bi-directional CDW state without the d-wave form factor^. 
The sign of the Nernst coefficient, on the other hand, is not di¬ 
rectly determined by the sign of the dominant charge carriers 
and thus may or may not be the same as the sign of the Hall 
and Seebeck coefficients. 

In this paper we ask if the quasiparticle Nernst coeffi¬ 


cient in the mean field BDW state does indeed show a drop 
with decreasing temperature, with v jT eventually becoming 
negative as T 0 as seen in experiments. We consider 
a two-dimensional (2D) bi-directional Q\ = (27r/3,0) and 
Q 2 = (0,27r/3) BDW state in mean field theory (valid for 
temperatures T < Tbdw and magnetic fields high enough 
to eliminate the superconductivity) and investigate the quasi¬ 
particle Nernst coefficient as functions of temperature and 
hole doping appropriate for the underdoped regime of the 
cuprates. Although the experimental evidence is that for a 
slight incommensuration in the BDW scattering vectors (i.e., 
qi ^ q 2 ^ 0.31) in this paper we work with a commensurate 
BDW for simplicity (i.e. we take qi = q 2 = 0.33, corre¬ 
sponding to charge modulations with periodicity of three lat¬ 
tice vectors). We find that, below the BDW transition temper¬ 
ature and in the appropriate regime of hole doping, the Fermi 
surface topology changes from a large hole-like Fermi sur¬ 
face at higher doping (where there is no BDW) to small Fermi 
surface pockets at lower doping. A similar fermi surface re¬ 
construction in terms of a CDW state was recently assumed 
to explain the low frequency o f qua ntum oscillations in the 
pseudogap phase of the cuprateJ^^^. We find that the quasi¬ 
particle Nernst coefficient in the mean field BDW state does 
indeed show a drop with decreasing temperature, with v jT 
eventually becoming negative as T ^ 0, as seen in experi¬ 
ments. 

This paper is organized as follows: In Sec. II, we consider 
the Hamiltonian for the BDW state and examine the energy 
spectrum and the reconstruction of the Fermi surface. In Sec. 
Ill, we define the quasiparticle transport coefficients which 
we compute numerically using Boltzmann semiclassical equa¬ 
tions. Sec. IV and V are devoted to analytical calculations of 
the Nernst coefficient in the limit of small order parameter and 
magnetic fields and the breakdown of Sondheimer cancella¬ 
tion. In Sec VI, we present our numerical results for Seebeck, 
Hall and the Nernst coefficient for the BDW state and show 
that they all become negative at low temperatures. We end 
with summary and conclusion in Sec. VII. Some analytic ex¬ 
pressions and formulas have been relegated to the appendix. 


II. MODEL AND FORMALISM 

In a mean-field picture, the Hamiltonian describing a den¬ 
sity wave ordered state can be written as, 

= Y, [W^(k)4+Q + /I.C.], (1) 

k,Q,cr 

where W{k.) is the order parameter which can in general de¬ 
scribe a charge, orbital current, or a bond density wave in 
cuprates depending on the form factor IV(k). The operator 
cr creates an electron of spin a with momentum k, and Q 
denotes the modulation wave vector. Charge modulations with 
a periodicity of 1/5 {1/5 integer) lattice vectors describes a 
commensurate CDW state. The modulations can be given by a 
uni-directional modulation Qi = 27r((5, 0) or Q 2 = 27r(0, (5) 
or a superimposition of the two wave-vectors in which case 







3 


♦ 4^4 ^9 41 


xlO 

1^ 


■ --2 


x/a 


FIG. 1. (color online) Bond order modulation in real space over a 
spacing of three lattice vectors. Plot of I'lpBDwix, y)\^ — y)\^ 

in real space of lattice constant a, where i^{x,y)I^BDw{x,y) is 
the lowest energy wave function without/with bond order parameter 
W = 0 . 2 eV. 


the CDW is bi-directional. A modulation wave-vector of type 
Q = 27r{6,6) can describe a third variant of the same CDW 
state. The functional dependence of the form factor W (k) and 
the modulation vector Q distinguish different density wave 
states for example Q = (tt, tt) and W (k) = (cos — cos ky) 
is the well known staggered flux or tZ-density wave (DDW) 
statJ^. 

The tight binding description for electrons on a two- 
dimensional square lattice of unit lattice constant is given by 
the energy dispersion relation 

Ck = —2ti (cos kx + cos ky) + 4^2 cos kx cos ky 

— 2ts (cos 2kx + cos 2ky)^ (2) 

where ti, t 2 and ts are the nearest neighbor, next-nearest 
neighbor and next-to-next-neighbor hopping parameters. For 
all numerical calculations, we chose the parameters ti = 1.0 


eV, t 2 = 0.3ti and ts = 0.1^2 which reproduce the non¬ 
interacting Fermi surface (see Fig. [^. We now focus on the 
Hamiltoniain for the BDW state. The following real space 
mean field Hamiltonian couples fermions to the bond ordeJSI 

HbDW = E [Wa (e*Q-('-+-/ 2 ) + eiQ-(rW 2 )) 

r,a,cr 

+ h.c], (3) 

where in the sum r denotes the lattice sites, the vector a rep¬ 
resents all the nearest neighbors vectors. The operator Cr,cr 
annihilates an electron of spin a at the site r. Recent ex¬ 
periments suggest that the bond order rese mbles the d- 
wave form factor W±x = —W±y = The vectors 

Qi = 27r((5, 0) and Q 2 = 27r(0, (5) describe the periodic mod¬ 
ulation of the bond order where ^ = 1 /3 indicating a com¬ 
mensurate BDW order with periodicity of three lattice vectors. 
The bond density wave order effectively redefines the hopping 
amplitude ti modulating it spatially. Fig. shows the bond 
order modulation in real space for the chosen modulation vec¬ 
tors. It is useful to Fourier transform Eq. ^ and rewrite the 
equation in momentum space: 

HbDw{^) = Wo E[(cOsfc3; - COsfcj^)4_^Q^/2,<TCk-Qi/2,<T 
k,cr 

4+Q2/2,<T^k-Q2/2,<T] + h.C, (4) 

where Ck,cr is the annihilation operator for an electron of mo¬ 
mentum k and spin a. The total Hamiltonian Hmf for the 
system is 


Hmf = ^ ekcj^ ^Ck,cr + Hbdw, (5) 

k,cr 

which can be expressed in terms of a nine component operator 
^'k,<T as 


Hmf = E ^k,a^(k)^'k,a (6) 

\z.^RBZ,a 

where RBZ is the reduced Brillouin zone (—tt/S < kx < 
7r/3, —7r/3 < ky < f/3) and i^(k) is 
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Details of the non-zero entries Wij in the Hamiltonian matrix trix imposes the condition that Wji = w^j. Diagonalizing the 
are given in the appendix. Note that hermiticity of the ma- Hamiltonian H(k) in Eq. we obtain the energy eigenval- 
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ues En(k) and the corresponding eigenvectors. Fig. shows 
the relevant bands of Hamiltonian Hmf near the chemical 
potential out of a total of 9 bands. One notes the presence of 
an electron pocket centered at (tt/S, 7r/3) and a hole pocket at 
(tt/S, 0) and symmetry related points, which are also depicted 
in the reconstructed Fermi surface in Fig. 

In Fig. 1^ we plot the electron spectral function for the 
BDW Hamiltonian. The electron spectral function A{uj,'k) 
is given by 

-4(w,k) =-IlmGi-et(w,k), (8) 

TT 

where Gi-et(<^ 5 k) is the retarded Green’s function for the 
Hamiltonian. A{uj,'k) essentially maps out the Fermi surface 
as it should be observed in ARPES experiments. In contrast 
to the Fermi surface plot in Fig. the electron spectral func¬ 
tion is not 27r/3 periodic in kx and ky, but it is weighted by 
the coherence factors at each point on the Brillouin zone^. 
A very similar ARPES spectral function for the BDW phase 
with slight incommensuration has recently appeared in Ref. 

El. 

The hole doping in the cuprates is conventionally counted 
from half-filling, i.e., one electron per Cu atom. If n denotes 
the fraction of occupied number of states in the Brillouin zone 
then the doping p = I — 2n. The fraction n is calculated as 

n= ^ /(^n(k)), (9) 

n,\i.ERBZ 

where/(F^n(k)) = 1/(1is the Fermi distribu¬ 
tion function which at zero temperature is simply a step func¬ 
tion 0(/i — F^^(k)). We find that p behaves linearly with dop¬ 
ing p. Half filling (p = 0) is evaluated to be at /i = —0.7055ti 
and a doping ofp = 12.5% is found at /i = —l.OOlfiti. 


III. QUASIPARTICLE TRANSPORT COEFFICIENTS 

The sign of the transport coefficients such as the Hall 
and Seebeck coefficients reveals information about the carrier 
types (electron or holes) and also the underlying Fermi sur¬ 
face. We use the formalism of linear response theory to cal¬ 
culate the Hall, Seebeck and Nernst coefficients for the BDW 
state. The charge current J and the thermal current Q can be 
related to the electric field E and the temperature gradient VT 
as 



The use of three conductivity tensors d, k and a is sufficient 
to relate thermal and electrical effects. The diagonal compo¬ 
nents of the matrix in Eq[^give us the electrical and thermal 
conductivity while a interrelates the thermal current and the 
charge current to electric field and the temperature gradient 
respectively. Applying a temperature gradient VT across the 
X axis of the sample, an electric field Ex is generated across 
given by Ex = axx'^xT and the Seebeck coefficient S 
is defined as S' = ocxxl^^xx^ Applying a magnetic field in the 



FIG. 2. (color online) Energy bands of BDW state for W = 0.22ey 
(a) along the path (0,0) to (2/3, 0), (b) along (2/3, 0) to (2/3, 2/3), 
(c) from (2/3, 2/3) to (0,0) (d) Energy bands in the limit of IF ^0. 
The solid black line indicates the chemical potential corresponding 
to a doping of 11%. 



FIG. 3. (color online) Fermi surface reconstruction for BDW state 
for p = 11% consisting of an electron pocket at (7r/3,7r/3) and hole 
pockets at (7r/3,0) and (0,7r/3). The green contour shows the Fermi 
surface without the BDW order parameter i.e. IF = 0. The smallest 
square at the center enclosed by dashed line is the reduced Brillouin 
zone (RBZ) appropriate for the BDW state. 


perpendicular direction now generates a Hall current Jy and 
the Hall coefficient is given by Rh = cFxy!cFxx^yy 

The Nernst effect measures transverse electrical response 
to a thermal gradient in the absence of a charge current i.e. 
Ey = dTIdx, where ^ is the Nernst coefficient and we 

apply —dTjdx thermal gradient along the x direction which 
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FIG. 4. Electron spectral function A(uj = 0, k) in the presence of 
bi-directional bond order Qi = (27r/3,0) and Q 2 = (0,27r/3) at 
a doping value of p = 10%. Fermi surface reconstruction due to 
BDW results in the formation of small electron and hole-like pockets 
which are also observed in Fig. The electron spectral function, 
unlike the bare Fermi surface, is weighted at each point of the BZ by 
coherence factors and therefore in general is not a periodic function 
of Qi or Q 2 . Similar spectral function for the BDW phase has also 
been reported earlieP^. 


is the appropriate experimental convention. From Eq. p^it 
follows that the Nemst coefficient is 

Q _ Ey _ ^xy^xx ^xx^xy 

{-dT/dx) <^xx + '^xv 

For magnetic field B pointing in the 2 ; direction, we redefine 
the Nernst coefficient to be = i9/B. The quantity vjT 
is the one which is determined experimentally. It is impor¬ 
tant to clarify the sign convention of the Nernst coefficient 
chosen here according to which the sign of the supercon¬ 
ducting Nernst signal is opposite to that of standard textbook 
conventioJ3l. According to this convention, the sign of the 
Nernst signal of the vortices is positive when there is a nega¬ 
tive temperature gradient along the x axis. 

We employ the semi-classical Boltzmann equations ap¬ 
proach for the calculation of conductivities in the relaxation 
time approximatioiP^ with the bi-directional BDW modula¬ 
tion. 

«- = ^E/(12) 

Oixy = 

n 

[vly^y - v^vyy)]d^k ( 13 ) 

= , 14 ) 


n 

(viv^y - vivy^y)\d?\<i, (i5) 

where n is the band index, is the semi-classical quasi¬ 
particle velocity and . The integra¬ 

tion is restricted to RBZ and the energy eigenvalues E^n(k) 
are measured relative to the chemical potential. The fac¬ 
tor of 2 present in the numerators takes into account spin- 
degeneracy of the energy bands and r(k) is the scattering 
time which takes in to account interactions between quasi¬ 
particles and impurities, phonons and other quasiparticles. 
We point out that r(k) is assumed to be independent of en¬ 
ergy but we retain a possible momentum dependence which 
yields a positive Seebeck coefficient in the normal state con¬ 
sistent with experiments. We assumed the scattering time 
r(k) = (1 -f a{coskx + cos/c^))^, where the parameter a 
is chosen to be 0.4^. The precise functional form of r(/c) is 
unimportant, however, and any other momentum dependence 
of the scattering time that produces a positive sign of the See¬ 
beck coefficient at high temperatures works just as well. Note 
that an assumption of a momentum independent r results in 
a negative Seebeck coefficient in the normal statd^^^ (i.e., 
above the BDW transition temperature), inconsistent with 
experimentSo although the Nemst coefficient u/T is 
robust and negative in the BDW phase as T ^ 0 even with a 
momentum independent r, which is our central result in this 
paper, we retain a momentum dependent scattering time only 
to be consistent with the sign of the high temperature See¬ 
beck coefficient (which is not the focus of this workj^. The 
temperature dependence of the conductivities arises from the 
factor of derivative of the Fermi function df{E{'k))/dE{'k) 
which takes the form of a Dirac-Delta function at absolute 
zero. 

Note that in our results we have set 5 = r(k = 

(7r/2,7r/2)) = 1. By examining the definitions of a and a, 
we note that we can make r(k) and B dimensionless by re¬ 
placing r(k) ^ r(k)/ro and B —> B{eT^t\o? jk?). Here 
To = r(k = (7r/2,7r/2)) is a representative scattering time 
that has been set to 1 and for our calculations a ^ 3.9A 
is also taken to be 1. Choosing B = 1 io correspond to a 
physical B ^ 2T (higher values of B does not qualitatively 
change our results), we obtain that a mean-scattering time of 
To = 1 chosen in this work corresponds to a mean scattering 
rate ^ lOA". 


IV. ANALYSIS OF NERNST EFFECT IN BDW STATE 
USING HOT-SPOT MODEL 


Before we compute the low temperature Nemst coefficient 
numerically using Eqs. 1^ to JA let us first understand the 
contribution to Nernst effect from mean field BDW state us¬ 
ing the so-called “hot spot” model shown in Fig.|^ For weak 
BDW amplitude W, the BDW can only affect electrons at the 
Fermi surface by scattering by wave-vector Qi. The strongest 
effect of the BDW is felt at momenta k, where both the start¬ 
ing wave-vector k and the ending wave-vector k Qi are on 
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the Fermi-surface. The wave-vectors k on the Fermi-surface 
which satisfy this condition are referred to as hotspots. 

The finite temperature thermoelectric coefficients and con¬ 
ductivities can be written in terms of the zero-temperature 
conductances as 

± 

cTj=°(Ai) = 26^^/ dkr^iv^fSiE^ - ^i), (16) 


y9v^ 


dv^ 


S{Ek - li) 



In the limit where crjj~^{p) varies slowly on the scale of 
T such that dpcrj^""^ » the thermoelectric coeffi¬ 

cients are given by the Mott relation 


aij 


TT^ T daJj-° 
3 e 5/i 


( 21 ) 


In this low temperature limit, the expression for Nernst signal 
can be simplified as: 


_ TT^ T OQh 
3 e dp 


( 22 ) 


where &h = -^o is the Hall angle. Since the Hall angle 
&H changes on the scale of the Fermi energy, which is much 
larger than the temperature, the typical contribution of a metal 
to the Nernst coefficient is small. This is referred to as the 
Sondheimer cancellatioiJ^. 

However, in the presence of a BDW order singular con¬ 
tributions to can lead to sharp changes in the Hall 

angle Qh, that result in enhancement of the Nernst coeffi¬ 
cient. To compute the effects of the BDW with wave-vector 
Q on T = 0 conductivity (i.e. Eq. [T^ [T^, it is convenient to 


change the integration variables from k^^ky to U = 


mdV = 


ek —ek+Q 


. The resulting conductivity is written as: 


= E / 

X 5{E - p) (23) 

E{U, V) = U + -f (a = ±1) (24) 

where (detailed expression are given in the appendix) 

are band velocities that are related to derivatives of the energy 

E. 
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FIG. 5. (color online) The ‘hot spot’ is defined as the point where 
one piece of unperturbed Fermi surface intersect with another piece 
when it is translated by the BDW wave-vector Qi. In this figure, blue 
lines are unperturbed Fermi surfaces, green ones are Fermi surfaces 
translated by vector Qi, the intersection between blue and green line 
at the tip of the arrow is the ‘hot spot’. 


To understand the role of the BDW order parameter W, we 
focus on the limit of a small order parameter W C \p\. The 
contributions to the thermoelectric and conductivity response 
are dominated by “hot spots” where the modification of Fermi 
surface by the BDW order is more dominant. In the small W 
limit, deviations from the Sondheimer relations can be under¬ 


stood in terms of singular terms in the integrand of Eq. [TC 
which are written as: 


dE _ V 
dV ~ “ 
d^E _ 

dV^ ~ '^( 1/2 _|_ ^ 2)1 ■ 


(25) 

(26) 


These terms develop singularity around the ‘hotspots’ where 
y ^ 0 in the presence of a small BDW order parameter W. 
The singularities at the hot-spot lead to linear in W contribu¬ 
tions to Eq. Since in this section we are interested in only 
linear W contribution to Eq.[^ the other terms in Eq.p^can 
be Taylor expanded at hot spots. 

Defining the BDW induced correction to as 


r / X W,T=0^ X W=0,T=0^ X 


(27) 


where ^ (p) is the conductivity when there is BDW in¬ 

duced gap, while cr]J“^’^“^(p) is the bare quantity. The 
leading order for <^crj“^(p) would be of 0{W), since results 
of 0{W^) is subtracted by the unperturbed ones. The non¬ 
singular contributions result in 0{W^) contributions which 
we will ignore. After simplifying Eq. (details in the ap¬ 
pendix), the leading order terms of da^yip) is written as 


^^xy{p) — —27r(-Fi -\r E 2 Gi -|- G 2 )W 

E7t{Es -E 4 EG 3 - G 4 )W (28) 


where the expressions for E and G are involved and given in 
the appendix. Similarly, we could obtain 5(7xx{pY 

ScTxxip) = -27r/iVF, (29) 
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FIG. 6. The change of electrical conductivities as a function of chem¬ 
ical potential /x. In this figure — 1 < /x < —0.8 is the range of 
chemical potential corresponding to the underdoping. 


where as before the expression for fi is given in the appendix. 
Expanding the conductivities for small W according to 
Eq.j^in Eq.j^we can write the linearized correction to the 
Nernst coefficient as 


50xy = 


3 e Ofi 


Ef 


tt^T d 
3 e dfji 




oO 


t 02 


Ef 


(30) 

(31) 


Using the parameters we set in the beginning of this paper, the 
numerical results for and da^y in the range of chemical 
potential to our interest are shown in EIG. Erom the plots, 
we see the dominant contribution to the Nernst coefficient in 
the chemical potential range —1 < /x < —0.7 comes from 
'§jl{^^xy) which leads to a positive Nernst signal. However, 
changes in parameters where the da^y curve shifts to lower 
chemical potentials will alter this sign to the experimentally 
consistent negative value. Such a shift can potentially be ob¬ 
tained by larger values of W beyond linear response, which 
produces a negative Nernst coefficient in a larger parameter 
range as we find in our numerical results described later. An¬ 
other interesting feature of Eig. is the divergence around 
/X ^ 0.3. This leads to a dramatic enhancement of the Nernst 
coefficient, which will be the subject of the next section. 


Nernst signal and a significant breakdown of the Sondheimer 
cancellation for specific structures of the fermi-surface. In 
Eig.|^ for a critical value of the chemical potential /Xc ^ 0.3 
both 6axy^ 5aXX appear to diverge. Examining Eq.[24| the di¬ 
vergence can be viewed as a result of the divergence in the 


Jacobian 


E{kx ) 
D{uy) 


_ D{kx^ ky) __2_ 

" D{U,V) ~ |v(k + Q)xv(k)|’ ^ 

which appears in Eq. 

To derive the form of the divergence as a function of chem¬ 
ical potential /x, we notice that for the Q = (Q, 0) symmetry 
of the BDW, 


x;^(ko + Q) = -'^a:(ko) Vy{ko + Q) = Vy{ko) (33) 

where ko is the position of hot spot. Extracting the divergent 
terms in 6aij we obtain: 


5crxy = ^ecBWqt 

5(Jxx = -^eWot 


'^a;,ko+Q 


'^cc,ko+Q 


Substituting the above, we obtain: 

T Va;,ko+Q 


= Wo 


3 

X 


7^ 

XX 






yy 


(34) 

(35) 


(36) 

(37) 


where Vy^x^^ = y 2Myy^(/x —/Xc), and Myy is the effec¬ 
tive mass at the hot spot. We see that there is a significant 
enhancement of Nernst signal around the critical chemical 
potential(/Xc). But the sign of the Nernst signal depends both 
on the bare Hall angle and the property of effective mass at 
hot spot and at critical chemical potential. 

This mechanism produces a large Nernst signal by break¬ 
ing down the Sondheimer cancellation when /x = /Xc where 
Eq. diverges. This occurs when the BDW wave-vector 
matches a nesting vector of the bare fermi surfaces. While 
the preceeding calculation is valid for small W and explains 
the breakdown of Sondheimer cancellation at a chemical po¬ 
tential beyond the usually accepted range in the underdoped 
regime, it is conceivable that larger values of W extend the 
range of breakdown of this cancellation. Besides, it has been 
proposecPS that strong correlation effects may lead to a low 
temperature electron pocket structure where the nesting would 
be similar to the required condition here i.e. /x ^ /Xc. Our 
result suggests that in addition to providing a natural expla¬ 
nation for the BDW wave instability with the experimentally 
observed wave-vector direction, this correlated statJSl would 
provide a mechanism for breakdown of Sondheimer cancella¬ 
tion even at small W. 


V. LARGE NERNST SIGNAL FOR SMALL W 

An interesting feature of the linear response results is a di¬ 
vergence in the conductance shifts 5a m Eig. Since the 
Nernst signal is proportional to the derivative of Hall an¬ 
gle over chemical potential, the divergence suggests a giant 


VI. NUMERICAL RESULTS FOR NERNST COEFFICIENT 
IN BDW STATE 

Eor numerical evaluation of the transport coefficients we 
choose values of the various parameters appropriate in the 
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T(K) 


FIG. 7. (color online) Plot of transport coefficients vs. temperature 
in arbitrary units for three different doping values showing an en¬ 
hanced negative signal in the underdoped regime (a) Seebeck coeffi¬ 
cient S/T, (b) Hall coefficient Rh and (c) Nemst coefficient u/T. 
The negative signals for doping values of p = 10% and p = 12.5% 
is ascribed to emergence of electron pockets due to Fermi surface 
reconstruction by BDW state. The Seebeck and Hall signal is posi¬ 
tive for a doping of p = 7% when W = 0 while Nernst coefficient 
shows a small positive signal close to zero (when compared to the 
large negative signal for the other two doping values). 



0.08 0.09 0.1 0.11 0.12 


P 



0.08 0.09 0.1 0.11 0.12 


P 



0.08 0.09 0.1 0.11 0.12 


P 

FIG. 8. (color online) Phase plot of transport coefficients in arbitrary 
units in the p — T space showing an enhanced negative signal in 
the underdoped regime and low temperatures (a) Seebeck coefficient 
S/T, (b) Hall coefficient Rh and (c) Nernst coefficient u/T. The 
dark blue area indicates the region of maximum negative response. 


pseudogap phase. It is important to emphasize that our re¬ 
sults (specifically the qualitative temperature and doping de¬ 
pendencies of the transport coefficients such as the Nernst co¬ 
efficient) are completely robust against variations of the nu¬ 
merical values and functional forms of the various parame¬ 
ters. The functional dependence of the BDW order param¬ 
eter on T,p etc is chosen to qualitatively mimic the exper¬ 
imental trend, and is by no means meant to produce quan¬ 
titatively accurate results for the Nernst and other transport 
coefficients in the pseudogap phase. We chose the bond 
density order parameter W = 0.22 eV at zero tempera¬ 
ture and at doping p = 0.125 (12.5% or 1/8 hole doping) 
which is also set as the upper critical doping value Pup, for 
bond order to survive. For any finite value of doping below 
Pup = 0.125, we assume a mean field doping dependence 
of W{p) to be W\{p - piow)/{Pup - Piow)\^^‘^, where we 
chose Plow = 0.085 as the lower cutoff. For any other doping 
range W{p) = 0. The bond order parameter is also assumed 
to scale mean field like with temperature below Tbdw{p) as 


W (p, T) = W (p) a /|1 — T/Tbdw{p)\- From experimental 
fit, Tbdw{p) is chosen to be 142 (p — p/)^'^, which gives the 
critical temperature in Kelvins. For T > Tbdw{p)^ W{p,T) 
is again chosen to be zero. 

Eq. [T^ to [T^ were used to calculate normal state conduc¬ 


tivities (a^x, (^xy) at a give temperature T and 

doping value p. Thus the T dependence of the the Seebeck 
{S = OLxxl(yxx). Hall {Rh = cFxylcJxxCTyy) and Nemst co¬ 
efficients {P = {axyCFxx - OLxxCrxy)l{cFlx + ^ly)) is Calcu¬ 
lated. Fig. 1^ shows normalized Seebeck and Nernst coeffi¬ 
cients {S/T and y/T) and Hall coefficient {Rh) as a function 
of temperature for three different doping values. For p = 10% 
and p = 12.5%, we observe negative coefficients ascribed to 
electron like pockets due to BDW order but the signal remains 
positive for p = 7% {p < pi) doping when TF = 0. Fig. 
shows a phase space plot of Seebeck, Hall and Nernst coeffi¬ 
cients in the p — T phase space, where the region of enhanced 
negative response can be visualized to be in the pseudogap 
regime of low temperature and low doping. 


VII. CONCLUSIONS 

In this work we studied the normal state of high Tc cuprates 
i.e in the absence of superconductivity when a large magnetic 
field is applied. Starting with a mean field Hamiltonian for 
bi-directional BDW order with wave-vectors Qi = (27r/3,0) 
and Q 2 = (0, 27r/3), we observed reconstruction of the Fermi 
surface from being large hole-like at higher doping (when 
there is no BDW order) to the appearance of small electron¬ 
like and hole-like pockets in the doping regime appropriate for 
the BDW state, which results from breaking of lattice trans¬ 
lational symmetry. The normal state Nernst effect is impor¬ 
tant to understand the Fermi surface topology of cuprates in 
the underdoped regime. The enhancement and the negative 
sign of the low temperature Nernst signal experimentally ob¬ 
served in the pseudogap phase of cuprates is the main focus 
of this work. The Nernst response typically vanishes for con¬ 
ventional metals due to Sondheimer cancellation, but this can- 
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cellation breaks down in the presence of magnetic field at low 
temperatures due to the presence of a BDW order parame¬ 
ter. In addition to providing analytical understanding for the 
breakdown of the Sondheimer’s cancellation in the presence 
of the BDW order parameter using the hot spot model, we 
numerically calculated all the three thermoelectric transport 
coefficients, namely Hall, Seebeck and Nemst coefficients, in 
the semi-classical Boltzmann approximation. At the temper¬ 
ature scale T < Tbdw^ observed a negative Nernst coef¬ 
ficient in the underdoped regime. This low temperature neg¬ 
ative Nernst response is reminiscent of a similar response of 
other two transport coefficients, namely the Hall and Seebeck 
coefficients. Though the negative sign of the Hall and Seebeck 
coefficients can be ascribed to the appearance of electron-like 
pockets which appear on the Fermi surface, the sign of the 
Nernst coefficient is not directly determined by the sign of the 
dominant charge carriers and depends on the detailed Fermi 
surface topology. 
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While calculating the leading order of we substi¬ 

tute in the expressions between eqs. | A1 1 and [A 11 1 into Eq. 
and perform the integration, and keep only the linear order 
terms in W. The result is Eq.|^ where. 


16 


Flip, 0) = rV(C/, F) iU,VyV,y + V,UyV,y + V,VyU,y) | (p.Q) 

(A12) 

F2{p,0) = A (T^JU.UyV.Vy) |(p,0) (A13) 

F3{P,0) = A [T^J{U,VyV,Vy+V,UyV,Vy)] |(p,0) 

(A14) 

F^ip,0) = A [r^J{V,VyV,Vy)] |(,,0) (A15) 


Appendix A: Expansion of Nernst coefficient 


The explicit expression for band velocity Va is 
^ dE dU dE dV dE 


V = 


dkr 


= U, + VA^\ 


dk^dU dk^dV 

w) 


vy = U, + V,{^ 


The inverse of effective mass terms are: 
dv^ 
dkx 
.... dv^ 


= Ux 


■Vx: 


~ ^yy ^yy 

_ 

— ^xy + ^xy 


'dE\ 


(k: 

.29^ F 


) + 

r 

’ dV^ 

dE\ 


iVy) 


w) 

+ 

dV^ 

dE^ 



. d^E 

W ) 

1 + 

VxVygv 2 ^ 


(Al) 

(A2) 

(A3) 

(A4) 

(A5) 


where 


Here, we have assumed the integral limits of V to be 
V/W -> Aoc. The above expression only represent the first 
part in the square bracket of the first equation in Eq. Eor 
the second part, we just change the derivative variable from 
xyxy to xxyy, a minus in the front and thus obtain G’s. When 
we consider Q^, because of the C 4 rotational symmetry, we 
simply interchange x and y. 


Similarly, we get the expression for 6 axx 


in Eq. 16 where 


fi{p,0) = tJ{V,)X, 0 ) (A16) 


E 

Uxx 
Uxy 


u±^/W+w^ 

dU _ dU 
^~Wy 
dV dV 

dkx ^ dky 

dkx [dU dkx [dV J 

dky \dU dky \dV J 

dU fdUy\ dV fdUy\ 

Wy \~W J^Wy / 


(A6) 

(A7) 

(A8) 

(A9) 

(AlO) 

(All) 


Appendix B: Hamiltonian matrix elements 


The non zero elements of the 9 component Hamiltonian are 
specifically given by: 
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Wi 2 = Wq {cos{kx ^ tt/S) — cosky) 

wis = Wq (cos {kx — tt/S) — cos ky) 

1 C 14 = Wq {cos kx — COS {ky IT/3)) 

Wi 7 = Wo {cos kx — cos {ky — tt/3)) 

W 23 = Wo (cos {kx + tt) — cos ky) 

'^25 = Wo (cos {kx + 27r/3) — cos {ky + 7r/3)) 

W28 = Wo (cos {kx + 27r/3) — cos {ky — 7r/3)) 

'^^36 = Wo (cos {kx — 27r/3) — cos {ky + 7r/3)) 

'^^39 = 11^0 (cos {kx — 27r/3) — cos {ky — 7r/3)) 

^C 45 = Wo (cos {kx + 7r/3) — cos {ky + 27r/3)) 

^C 46 = Wo (cos (/Ca^ — 7r/3) — COS {ky + 27r/3)) 

1 C 47 = VFo (cos kx — cos {ky + tt)) 

'^^56 = Wo (cos {kx + tt) — cos {ky + 27r/3)) 
'^^58 = Wo (cos {kx + 27r/3) — cos{ky + tt)) 

Woo = Wo (cos {kx — 27r/3) — cos{ky + tt)) 
'^78 = Wo (cos (/Ca^ + '7r/3) — cos {ky — 27r/3)) 

W 70 = Wo (cos {kx — 7r/3) — cos {ky — 27r/3)) 

wso = Wo (cos {kx + tt) — cos {ky — 27r/3)) 
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